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FOREWORD 


This is a final report for the second phase of the program. Evaluation 
and Prediction of Long-Term Space Environmental Effects on Nonmetallic 
Material j, conducted by Martin Marietta for the National Aeronautics and 
Space Administration, Marshall Space Flight Center, under contract 
NAS8-33578. 

The program was conducted in the Mechanical Materials Engineering Sec- 
tion with Mohan Misra as Program Manager and Harold Papazian as Prin- 
cipal Investigator. Contributors to the program included W. J. Arbegast 
and C. E. Forsyth for material tests. 

Up to the time of his unfortunate demise, Charles Peacock of NASA-MSFC 
served as program monitor. Dr. Ray Gause and Ms. Ann Whitaker served as 
technical advisors. 
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INTRODUCTION 


Although nonmetallic materials are generally recognized as being sus- 
ceptible to the changes caused by the environments associated with long- 
life spacecraft, data concerning the effects of irradiation and long- 
term thermal-vacuum exposure on the specific mechanical, electrical, and 
thermal properties of such material have not been readily available. 
Because the performance and life of any structure, system or component 
directly depends on the integrity of the materials used in its con- 
struction, the ability of these materials to withstand the operating en- 
vironment is of utmost importance in meeting service life requirements. 


The degradation of any material or material s.y^tem, even if slow at 
normal -use temperatures, can lead to failure over a long period. For 
example, in a polymer system used as a structural member, the accumu- 
lated changes may eventually lead to loss of integrity. Real-time tests 
of materials that may be used in a system can be clearly impractical in 
time nd money. If short-time tests can be used to predict future be- 
havior, these tests can considerably enhance mission success. 


The primary emphasis of the preceding pftise of this program was on the 
effects of long-duration thermal-vacuum exposure and had its foundation 
in the Viking materials qualification program of some 10 years ago. The 
project established the requirements for materials suitable for the 
Viking spacecraft. Specimens of the same materials that had remained in 
a thermal-vacuum environment to the present time were tested to deter- 
mine whether any changes in material properties had occurred over the 
long-term exposure. In general the largest changes occurred in such 
electrical properties as dielectric strength, dielectric constant and 
volume resistivity. The mechanical properties showed only small 
changes — some increasing and others decreasing in strength. The pre- 
ceding phase thus established the important fact that thermal-vacuum ex- 
posure would not degrade the structural integrity of a spacecraft con- 
structed with nonmetallic material systems. 


The present phase of the program was primarily devoted to the effects 
of charged particle irradiation on the mechanical properties of gra- 
phite/epoxy composites. There is a paucity of literature dataw con- 
cerning the charged particle irradiation effects on such composites, and 
it is essential that this task be performed as soon as the accelerator 
becomes available. The test results obtained to date are the baselines 
to be used later for comparison with results after irradiation. In pre- 
paration for the future testing, Section V discusses in considerable de- 
tail the effects that may be observed. 


Considerable effort was devoted to development of models and accelerated 
test methods useful for predicting material behavior. Zhurkov's ap- 
proach to the strength of solids is a useful first step toward nuch pre- 
dictions. This is discussed in considerable detail showing the under- 
lying fundamentals of the approach as well as its application to ac- 
celerated testing of composites. 



PROGRAM OBJECTIVE 


The objective of this program is to continue the experimental ^valuation 
of changes in the functional properties of a number of nonmetallic mat- 
erials as a function of simulated space environments and to use such 
data to develop models for accelerated test methods useful for pre- 
dicting such behavioral changes. 


Another objective is to evaluate the effects of charged particle irradi 
ations on candidate space materials. 


III. TECHNICAL APPROACH 


A. TASKS 

Task 1 requires evaluation of the final 20 Viking specimens. These were 
to be irradiated and the results compared with the preirradiated results 
obtained in the first phase of the program. 


Task 2 requires testing graphite/epoxy specimens prepared by Hercules. 
Sixty-four in-plane shear specimens were to be tested— 20 baseline spec- 
imens, 20 specimens to be tested after exposure to thermal-vacuum, and 
24 specimens to be tested after irradiation of specimens that had under- 
gone the thermal-vacuum exposure. 


Task 3 requires testing of graphite/epoxy specimens prepared by MSFC 
that had undergone thermal-vacuum exposure for some six years at MSFC. 
They were of four different configurations for appropriate mechanical 
property tests. These were tensile, interlaminer short beam shear, 
flexure, and compression strength. Twenty-three specimens were to 
tested as a thermal-vacuum-exposed baseline and 41 specimens were to be 
tested after irradiation by MSFC. 


Task 4 requires testing uf thermal-vacuum-exposed specimens after 
charged particle irradiation. 


Task 5 requires comparison of pre- and postexposure data for the mat- 
erials previously described and an attempt tc develop analytical models 
for predicting the long-term environmental degradation of selected mat- 
erial properties. Expected irradiation effects and predictive testing 
are discussed in considerable detail. 


A total of 64 graphite/epoxy specimens were tested for baseline {preir- 
radiation) data, and a total of 84 specimens were to be tested after 
undergoing irradiation. Of the latter, 20 were Viking specimens and 64 
were to be graphite/epoxy for comparison with the preirradiation base- 
line results. 
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TEST PROCEDURES 


The most Important aspect of the MSFC program is the fact that once a 
specimen enters the thermal-vacuum environment, it is never again ex- 
posed to the ambient environment. All specimens are transferred be- 
tween the thermal-vacuum system, the irradiation chamber, and final 
transfer to the in situ vacuum chamber under a vacuum environment, ^ hi s 
is done by transferring canisters &ntaining specimens under vacuum to 
various systems and finally through the airlock of the test chamber. 

The importance of this approach is demonstrated in Section V.A. 


All tests were conducted according to appropriate ASTM test methods and 
are designated in the appropriate table headings in Section IV. 


C. SPECIMEN DESCRIPTIONS 


The Hercules graphite/epoxy specimens were ±45 degrees in-plane shear of 
4- and 8-ply AS4/3501-6 and HMS/35 01—6 materials. The specimens desig- 
nated in Table 2 (Section IV) as 004 and 005 were of AS4/3501-6 with 6k 
tows and 12-mil and 21-mil thickness, respectively. The specimens de- 
signated as 009 and 010 were of HMS/3501-6 with 12k tows and 12-mil and 
21-mil thickness, respectively. U) 


The graphite/epoxy composites, Type 37 etc, were prepared by MSFC some 
six years ago and remained under thermal-vacuum exposure through 
testing. At present the preparations for the various types are not 
available. Within the four configurations, the types selected for 
testing included those with enough specimens for meaningful comparisons 
between baseline and irradiation tests. 
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IV. RESULTS 


All values reported In the following tables are the averages of several 
specimens. The * values indicate the maximum spread from the average 
and are probably attributable to differences in specimen thicknesses, 
Only the average thickness is given in the tables. 


Table IV-1 shows the results for the Hercules graphite/epoxy specimens 
under In-plane shear. In all cases modulus increased after a year under 
thermal-vacuum exposure. (Five scmples of each specimen type were 
tested.) This is probably attributable to further curing and/or out- 
gassing of a pi ast i sizing agent of some kind. Although the change in 
ultimate strength of the AS4/3501-6 (004, 005) may not be significant, 
the changes for the HMS/3501-6 material (009, 010) are large enough to 
be considered significant. Table IV-2 shows the resin contents of these 
materials. The TGA method is much simpler than the standard acid de- 
composition method, and we consider it to be more nearly representative 
of the true resin content. Comparison of these results with those in 
Table IV— 1 shows that the decrease in ultimate for both materials is as- 
sociated with the smaller amount ef resin content. If this is real, a 
significant observation has been mude and suggests more testing in this 
direction for in-plane shear studies. 


Tables IV-3 through IV-6 presents results for the specimens prepared at 
MSFC that were under thermal-vacuum exposure for some six years. These 
results serve as a baseline for comparison with the postirradiation 
tests to be performed in the next phase of this work. 


Table IV-3 shows the results of tensile testing. Assuming the same 
resin was used and noting the difference in thickness, it must be 
concluded that two different fibers were used in the two types. 


Table IV— 4 shows the apparent shear strength of Type 37, and Table IV-5 
shows the flexure strength of Type 38. These will serve as baselines 
for later irradiation tests. 


Table IV-6 shows the results of compression testing. Some results with 
electron irradiation obtained during the last phase of the program have 
been included. Because epoxy resins undergo crosslinking under irra- 
diation, the behavior of Type 26 is as expected — an increase in strength 
after irradiation. Only two specimens of Type 38 were irradiated and 
the slight decrease in strength may not be significant. If much less 
resin were used than for Type 26, the decrease could also be attributed 
to that fact. 


In conclusion, a much larger number of specimens must be tested than has 
been possible to date, and much better quality control during specimen 
preparation and irradiation must be exercised. Similar comments have 
been made by others. 
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Table IV-1 Hercules Graphite/Epoxy Tensile Strength (ASTM D3410) 


Spsclmen 

Batalina 

1-Year Thermal-Vacuum Exposure 

% Change of Averages 


Ultimata, ksi 

Modulus x 10 6 

Ultimata, ksi 

Modulus x 10® 

Ultimata 

Modulus 

AS4/3501-6 







004 

4-ply/l2.3 mil* 

2.37 

2,49 

+3 

+5 

005 

8-ply/21.3 mils 


2.62 

24.3^ 

3,02 


+15 

HMS/3501-6 

' m *\a 


11 ‘ 5+ .?,6 




009 

4-ply/12 mils 

2.33 

3,00 

+11 

+29 

010 

0-ply/21.3 mils 

12 ’ 7 tl|3 

2.51 

"*r* 

2.55 

•9 

+2 


Table IV -2 Resin Content Table IV~3 MSEC Graphite/Epoxy 1 ensile 


Strength (ASTM D2990) 


Specimen 

Ultimate, ksi 

Modulus x 10® 

% Elongation 

Typo 36 
<62 mils) 

27 1 +9 ' C 

. 3,3 

6 04 + ®‘® 
-0.77 

o ni +0 ' 4 
0,81 >0.46 

Typo 41 
(36 mils) 

7+10.7 

36 ’ 7 -7.3 

6 14+0.39 
.0.63 

- 07 +0.63 
1,07 +3.33 


Series 

Harculas Data, % 

TGA (Martin Marietta), % 

004 

33.4 

27.4 

005 

26.9 

22.8 

009 

28.4 

28.0 

010 

24.0 

25.4 


Table IV-4 

Mb FC Graphite/Epoxy Interlaminar 
Short Beam Shear Strength (ASTM 
D2344) 6 


Specimen 

Apparent Shear Strength, ksi 

Typo 37 
(88 mils) 

B C7+0.45 
6,67 -0.23 


Table IV-5 

MSFC Graphite/Epoxy Flexure 
Strength (ASTM D790) 


Spoclmen 

Maximum Fiber Strength, ksi 

Typo 38 
(120 mils) 

kr +7.4 
56 >1.2 


Table IV-6 

MSFC Graphite/Epoxy Compression Strength (ASTM 
D3410) 



Compression Strength, ksi 

Specimen 

6-Year Thermal-Vacuum Exposure 

After Irradiation 

Typo 24 
(80 mils) 

33.9 tg 


Type 26 
(90 mils) 

29.7 

46, T^ 0 

Typo 37 
(90 mils) 

44.2 + . 1 7 °| 


Typo 38 
(100 mils) 

K7 0+7.4 
62>2 -3.4 

47 5 +2,8 
-2.8 

Typo 40 
(61 mils) 

. n r +12.4 
40 - 6 - 10.8 


Typo 42 
(62 mils) 
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V. 


ANALYSIS 


A. IRRADIATION DAMAGE ESTIMATES 


Although the experimental portion cf irradiation damage was not com- 
pleted, the expected results were assessed in preparation for the test 
results. 


The most important effDcts of irradiation in polymers are crosslinking 
and degradation. If crosslinking predominates, the ultimate effect of 
irradiation will be to produce a network polymer in which all molecules 
are joined to each other. If degradation predominates, the molecules 
become smaller and smaller as irradiation proceeds and the material 
loses its polymeric properties. 


The degree of unsaturation of a polymer can also change on irradiation. 
Polymers that are initially highly unsaturated tend to become less so on 
irradiation, while with initially saturated polymers the amount of un- 
saturation increases with dose. The formation of an unsaturated double 
bond will tend to labilize the bonds (free to move) in the allyl posi- 
tion, 1 ^iing to a conjugated double bond system. Such an effect ac- 
counts i w the coloration of some polymers on irradiation. Another im- 
port sat change is the formation of gas during irradiation. The amount 
of g^s produced depends on the polymer, dose, temperature, type of ir- 
radiation, etc. 


The free radicals (highly labile molecular fragments) formed during ir- 
radiation^ must play an important part in chemical changes. These 
free radicals can be both short-lived and long-lived depending on the 
polymer. The response of polymers to irradiation is markedly affected 
by the presence of other substances, most probably from reaction with 
the free radicals produced. Oxygen produces a marked effect, usually a 
decrease in the tendency to crosslink, i.e., oxygen reaction with the 
broken bonds (of long-lived radicals) rather than crosslinking. Thus 
any irradiation study that allows the sample to be exposed to the am- 
bient environment after irradiation must be suspect. (The experimental 
methods of MSFC overcome such concerns because the specimens are never 
exposed to ambient.) The amount of oxygen initially present in the 
material causes only a small effect. Large effects can occur if the 
sample is thin enough and the dose rate is low enough for oxygen to dif- 
fuse during the irradiation. 


* • 
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To produce marked changes in most polymeric materials requires roughly 
one chemical change per molecule. A dose of R rads (1 rad = 100 ergs/ 
gm) on a molecule of M molecular weight deposits an average of 


(R x mo ergs/gm)| 


M gm mo1e~l 


i5*02 x 10 23 molecules mole^Tx 1.602 x 10* 2 ergs eV _1 y 


- 1.04 x 10“*°(RM) eV/molecule. 


Denotifig the chemical changes per 100 eV as G, the dose needed to ac- 
hieve a change is 

R x (Ml/o.96 x 10101 * 0.96 x 10^ rads. 

\ G J\ M / 5M 


In all but exceptional cases G ranges between 1 and 5. Taking M « 2000 
and G = 1 yi^los R = 4.8 x 10 8 rads. For the worst case with G = 5, 

R x C„6 x 10' rads. In any long-duration flight these doses are 
easily obtained. 


Figure 1 shows the penetration depth as a function of dose in a 30-year 
geostationary orbit for an organic composite. (4) It is represented by 
an exposure of approximately 10 11 rads at 0.5-um (0.02-mil) depth to 
low-energy protons and approximately 5 x 10 9 rads at 50-um (0.2-mil) 
depth due to electrons. As expected, energy deposition decreases with 
depth in the material and shows that a 30-year flight doses on the order 
of 10° rads would cause changes as deep as 100 mils in relatively 
thick composites. 


os? s 
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Penetration Depth, mils 


Figure V-l Penetration Depth As a Function of Dose in Geostationary Orbit 





Figure 2 shows a decay rate analysis (developed from data in Ref 5) of 
free radicals in an epoxy resin (TGDDM-DDS) at room temperature after 
electron irradiation. Clearly long-lived radicals are present. As pre- 
viously discussed, if after irradiation specimens are exposed to the am- 
bient environment before testing, the results must be questioned. The 
MSFC approach avoids this concern by never exposing specimens; the fi- 
gure confirms the importance of our approach. 


Only one cursory reference to proton irradiation was revealed by a li- 
terature search(l). After 756 hours at 10“ 6 torr with 4 suns UV radi- 
ation and 10 minutes 1-MeV proton irradiation (dose not given), HT-S/ 
X-904 showed no change in flexure, but shear decreased by 14 percent. 
GY-70/X-904 showed a 10 percent decrease in flexure but only a 1.5 per- 
cent decrease in shear. It should be emphasized that these changes may 
be associated with the UV radiation because no tests with only UV were 
performed. Using 1-MeV electrons instead of protons, a 12 percent de- 
crease in flexure and a 19 percent decrease in shear was observed with 
HT-S/X-904, while GY-70/X-904 showed a 21 percent, increase in flexure 
and a 3 percent increase in shear. 


The flexure properties of miniature (22 x 13 x 0.51-mm) unidirectional 
specimens of T300/5208 composite were exposed to electron irradi- 
ations' 3 ) with doses up to 8 x 10 9 rads. The strength increased 
slightly with dose (approximately 4 percent), while the modulus was un- 
affected. Statistically, however, the limits for the 95 percent confi- 
dence led to uncertainty in the location of the mean. The authors sug- 
gested the possibility of large, perhaps significant, changes in flexure 
and emphasized the need for quality control in specimen preparation and 
testing. 


These comments demonstrate the paucity of literature data concerning the 
effects on composites as well as the uncertainty of available results 
and the importance of the MSFC program. 


A further consideration (generally not discussed) is the radiation of a 
mechanical property while under load such as might occur during reorien- 
tation or change of orbit. Figure 3 shows such an effect on the creep 
rate of polystyrene developed from data given in Reference 6. It may be 
of value to consider such phenomena in some further phase of the pro- 
gram. 
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B. 


ZHURKOV'S APPROACH TO THE STRENGTH OF SOLIDS FOR PREDICTIVE MODELING 


The mechanisms and the laws governing the failure of brittle solids and 
polymers have features in common and understanding the phenomena can 
lead to successful development of accelerated tests for such materials. 
For long-duration performance of ...aterials, the designer needs to have 
data on how long his choice of material can be expected to last in a 
given situation. Real-time tests are clearly impractical in both time 
and money, and the obvious solution to the problem is to use an accept- 
able accelerated test to predict lifetime. 


The time dependence of the strength of solids and. its connection with 
the mechanism of failure was studied by Zhurkov! 7 ) and his coworkers. 

His phenomenological approach to delayed failure (creep rupture) is also 
known as kinetic rate theory. In the Soviet Union the Zhurkov method is 
unquestioned! 8 ) while in the West it has been suspect leading to ques- 
tioning of the method with such comments as "... a good bit of discre- 
tion is often required to force the data to intersect at a pole."! 9 ) 

The purpose of the present discussion is to minimize the discretion. 

The Zhurkov method has modifications under certain conditions that have 
not been recognized. In some materials the common pole 0 , which is 
obtained from log tvs a data at several temperatures, is much larger 
than t 0 ~ 10~ 13 seconds obtained with materials that conform to the 
usual Zhurkov method. Zhurkov associated r 0 with the vibration fre- 
quency of the bonded atoms in the solid. In other instances experi- 
mental data do not result in a common pole at all. The following dis- 
cussion shows how these materials can be brought into the Zhurkov frame- 
work. 

From an applied viewpoint, the primary value of the Zhurkov method is 
that it may be used to devise accelerated testing of materials. Here 
such tests are outlined for creep rupture within the usual method and 
its modifications. The following section shows that accelerated tests 
for composites under static, dynamic and cyclic loads are also possible. 

Zhurkov developed his creep rupture method by successfully applying his 
relationships to a variety of materials (metals and nonmetals) and their 
time dependence of strength under unidirectional tension. 

The relationships and definitions of Zhurkov are: 


[1] f = Ae~“° (at constant temperature) 

[2] r = r oe U*/KT (at fixed a) 
where 

r « time to rupture, 
a ss constant, 

A = constant, 

a *» applied uniaxial stress, 
r 0 = constant, 

U* = U 0 - y 0 (energy), 

U 0 = an activation energy, 

Y = constant, 
k = Boltzmann's constant, 

T = absolute temperature. 

When a material behaves according to Equations [1] and [2], js con- 
sidered to be the period of natural oscillation, i.e., t 0 = 1/v where 
v is the vibrational frequency of the bonded atoms; U 0 is the acti- 
vation energy of the elementary act of failure in the absence of stress; 
and y is a coefficient that depends on the structure of the material. 
Equations [1] and [2] are related through 


[3a] a - y/kT 

[3b] A = r 0 exp(U 0 /kT). 


Figure 4 shows a schematic of results obtained by applying Equation [1] 
to creep rupture data. Such results are typical of those obtained by 
Zhurkov for a variety of materials. From such results, the constants A 
and a are obtained at each temperature that is studied. 


Figure 5 shows a schematic of results obtained by applying Equation [2] 
to such data where the slopes equal U*/2.303k. A U* is obtained for 
each a. Then a plot of U* vs a yields U* = - ya + U 0 , where the slope 
is y and the intercept Uq. It should be noted that y may also be ob- 
tained from Equation [3a j- 
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Data similar to Figures 4 and 5 for viscose fiber may be found in 
Chapter 1 of Reference 8 and are analyzed in complete accord with the 
Zhurkov method. 


Thus for an accelerated test of creep rupture, at some temperature of 
interest three or four points at relatively high stresses (short rupture 
times) will permit predictions of results at much lower stresses, i.e., 
at much longer times to rapture. Such data taken at two or three easily 
accessible temperatures should predict results for desired stress at, 
say, very low temperatures. Such applications are the claimed useful- 
ness of the method, and have indeed been applied to many materials. (w 


In Chapter 1 of Reference 8 the time dependence of the strength of Cell- 
uloid at four temperatures (20, 70, 85 and 100°C) is presented (without 
discussion) in a figure similar to Figure 4 here. However, the pole 
r Qt occurs near 5 x 10"^ seconds — some nine orders of magnitude re- 
moved from normal vibrational frequencies for well-behaved Zhurkov ma- 
terials. From these data then, Equation [3b] is A = 5 x 10 ~ 4 exp (U /RT). 
Here and below, R will be used rather than k and the energy terms will 
be in calories. 


The data for Celluloid were replotted as in Figure 5. From the slopes 
the U* were obtained as a function of a. As before this yields 
U* = -yer + U . For Celluloid the results are U* = -5.9 a + 62 (kcals). 
Application of Equation [3a] using the a's obtained from the original 
data give values of y only some 20 to 40 percent of 5.9. Clearly 
Celluloid does not behave "normally" for the Zhurkov method where the 
y's from Equation [3a] or [3b] yield the same value. 


A somewhat different analysis of the Celluloid data reveals thermo- 
dynamic functions not apparent in the preceding discussion. From Equa- 
tion [1], for any fixed , the constant A may be obtained as a function 
of temperature. The results are presented in Figure 6 from which A = 
5.6 x 10~ 3 ^ exp (63400/RT). The preexponential factor is 30 orders of 
magnitude from that obtained above (i.e., 5 x 10“ 4 ). With Zhurkov 
materials, the preexponential factors do not differ for a similar anal- 
ysis. The energy terms (U 0 ~ 62 ~ 63.4 kcals) are in good agreement 
considering the fact that data were extracted from the rather small fi- 
gure of Reference 8. 


Using 

these 


the value i: n =5 x 10“ 4 and A's from Figure 6 and inserting 

. I r* J ^ n I n it y I I / t m i . i * * . 


distinguish from the above-determined U 
function of temperature shown in 
cals. It appears then that Uq = 
here that this is of the form aF 


Figure 7. 
-130T + 

<= aH - 


into Equation [3b], A = r Q exp (Uq/RT), where the prime is* to 

U A may be obtained as a 
, where = -130T + 62000 
U 0 . It should be noted 
TaS. 
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For these two different A's, i.e,, A « 5 x 10- 4 exp (U 0 /RT) and A « 

5,6 x lQ-^exp (Ug/RT), Zhurkov's Equation [1] yields two equations 
that must be resolved: 

t a 5 x KT^eUg/RTg-acr 
and 

t* 5.6 X lQ-34 e U 0 /RT e -atf 
Substituting for Ug * -130T + U 0 
ln^a -130/R + Ug/RT - aa + lntg 
where ^g should be the 5xl0~ 4 value, and 
In Tn Ug/RT - aa + ln(5.6 X 10~ 34 ). 

Equating and solving for *tg 

x o “tge^O/R « 5.6 x 10“ 34 (3.27 x 10 28 ) » 1.83 x 10~ 5 , 

which is in reasonable agreement with 5 x 10“ 4 considering the avail- 
able data. 


It should be noted that the exp (130/R) Is an entropy of activa- 
tion! 10 ) . This also follows from the analogy between U 0 and aF men- 
tioned above. Then U 0 is the enthalpy of activation aH and 

i* i* 

t=5 x io-4 e AS/R e ^H/RT e - aa 
or 

t 

t „ T e *F /RTe -a a 
t = TQeC^F^-ya/RTj^ 

which is of course Zhurkov's Equation [c.']„ Thus for "normal" Zhurkov 
materials, where Up is obtained by extrapolation of U(cr) to a = 0, 

U 0 is to be Identified with the free energy of a standard state. 

In air, Lucite! 11 ) and Delrin! 12 ) behave in a manner similar to 
Celluloid, i.e., the pole depicted in Figure 4 is also some 9 to 10 
orders of magnitude removed from Zhurkov-type materials( Tn ~ 10~ 13 ) 
They can be assessed in the same manner as above for Celluloid. 


Accelerated tests for these kinds of materials can be carried out in the 
same manner as noted above for "normal" Zhurkov materials. 


When Lucite and Delrin are tested under vacuum, the data for log vs a 
do not intersect at a common 0 as depicted in Figure 4. A series of 
parallel lines for the various temperature are observed. Similar re- 
sults are observed for cellulose, acetate and polystyrene for both air 
and vacuum tests( 11 ). 


Figure 5 then is modified to Figure 8 when log is considered as a 
function of temperature (Zhurkov 's Equation [2]). 


The existence of parallel slopes for log 
considering Figure 8: 
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The following analysis shows that y * 0. However, in filled polymers 
such as graphite/epoxy composites, it is shown in the next section shows 
that y a f(c). 


Because log r vs a have parallel slopes, a is a constant, and for Lucite 
o » 3.36 at any temperature. From results similar to those depicted in 
Figure 8, U* a 42.8 kcals. An analysis similar to that depicted in Fi- 
gure 5 yields 

A = 8.2 x 10-2° e u */ RT . 


An entropy term on the order of e* 3,8 ^ will bring the preexponential 
factor in line with a Zhurkov r 0 . 

Analysis for r 0 vs a (shown in Fig. 9) yields 

r 0 = 8.2 x 10“2Q e“ 3,3 ^ 0 . 


Because a = 3.36, then 
Tq p 8.2 x IQ" 20 e~ acr . 


Inserting A into Equation [1] and into [2] and equating their 
logarithms shows U* = U 0 - ya 9 and because U* is constant (Fig. 8) it 
must be concluded that y = 0 for materials showing parallel slopes for 
log rvs 0 at various temperatures. 


Finally, for results showing parallel slopes, Zhurkov's Equation [2] is 
simply 

r = r 0 eV R T. 


This conclusion was reached earlier^ 11 ) without detailed analysis. 
It has been shown^ 3 ) that for such materials U 0 is essentially 
identical to the activation energies of thermal decomposition. 


For these kinds of material behavior (i.e., parallel slopes), the ac- 
celerated tests are simpler to formulate because the temperature de- 
pendence can be obtained with fewer test points. 


Zhurkov claimed that his equations were verified for both metals and 
nonmetals. For nonmetals, at least, their modifications must be con- 
sidered when undertaking data analysis by his method. However, the ap- 
proach can be used to devise quite simple and time-saving accelerated 
tests and warrants more attention than it has received. 
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The discussion most clearly indicates the kinetic rate theory basis of 
the approach. Zhurkov's energy term U 0 is to be associated with a 
standard state where this state is for a « 0. This appears to be con- 
firmed by the association of U 0 with the activation energy of thermal 
decompositions' 11 /, where obviously a « 0. (It would be instructive 
to study thermal decompositions of materials under load.) Since the 
Arrhenius rate constant is given by kf = Z exp (-aE/RT), the constant 
A in Zhurkov's Equation [1] appears to be A - 1/ky with r 0 = 1/Z. 

From absolute rate theory the rate constant is given by 

k r = (kT/h)e AS /Re“ a ^ W 

■\nd therefore Z can be equated with (kT/h)exp ( aS^/R) . The kT/h - v 
second" 1 is associated with vibrations in an "activated complex." 

Then t q - 1/v (as defined by Zhurkov) ranges between (1.6 - 
1.2)10~ 13 . These considerations imply the absence of an entropy term 
in "normal" Zhurkov materials as the analysis has shown. For materials 
where r 0 £ 10“ 13 , the entropy term appears to be contained within 

r o- 

Comparison of the behavior in air and in vacuum of Lucite^ 11 ) and 
Delrin^) shows that for some polymers y should be associated with 
surface phenomena as suggested in the next section for the graphite/ 
epoxy interface. If this conjecture is viable, then from studies of 
such materials it may be possible to determine the activation thermo- 
dynamics of crack growth. However the true role of y must await fur- 
ther work. 


C. ACCELERATED TESTING OF COMPOSITES 


Modern technology is making extensive use of composites, especially in 
aircraft and space structures. For long-duration performance, the de- 
signer needs to know how long a composite under load can be expected to 
last in a given situation. Real-time tests are clearly impractical in 
both time and money. The obvious solution to this problem is to use ac- 
celerated test methods to predict the lifetime of the composite. A re- 
view of the literature shows a paucity of work on accelerated testing 
in general and particularly for composites. Brinson and his 
coworkersl 13 ) have tried the time-temperature superposition method and 
Chiao' 14 ) has used the creep rupture method of Zhurkov. Both 
methods were developed before the advent of composite materials and have 
shown only modest success when applied to composites. The Zhurkov 
method's success had been modest because of a lack of understanding of 
its limitations and modifications. 


It is the purpose of this section to show that the Zhurkov method can 
also be applied to dynamic tension and to cyclic loading and provides a 
viable approach to accelerated testing of composites. 
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The following data (obtained from the literature) will show a straight- 
forward application of Zhurkov's Equation [1] for dynamic tension of 
glass fiber and cyclic loading for glass/polymer, metal matrix and grap- 
hite/epoxy composites. Temperature data were available for only the 
latter. 


Zhurkov developed his creep rupture method by successfully applying his 
relationships to over 50 different metals and polymers. The Russian 
literature indicates that his method is unquestioned there. Outside 
Russia the method has received little attention and, as developed in the 
following sections, we will show how unfortunate this oversight has 
been. 


The relationships of Zhurkov are as before 
T - Ae"“ a (at constant temperature) 
r- r 0 e u */ kT (at fixed a). 


According to Zhurkov when a material behaves according to Equations [1] 
and [2], r„ is considered to be the period of natural oscillation in 
the solid (~lcH-2 seconds) and y is a constant that depends on the 
structure of the material. Then Equations [1] and [2] are related 
through A = r 0 exp (U 0 /kT) and a = y/kT. 


For an accelerated test at some temperature of interest, a few data 
points at relatively high loads (i.e., short failure times) will permit 
evaluation of the constants in Equation [1]. Then Equation [1] can be 
used to predict failure times at low loads where the failure times are 
long. With high load tests at two, or at the most three, temperatures, 
the constants of Equation [2] can be evaluated. Failure times can then 
be predicted at any temperature and load. The results below show the 
efficacy of Equation [1] to a variety of composites. Temperature data 
were found only for the graphite/epoxy composites. 


Dynamic Tension 


Figure 10 shows results for Suprasil-2 glass fibers under static and 
dynamic tension. Zhurkov's Equation [1] may be applied to both con- 
ditions. Taking a few data points at relatively high loads (short rup- 
ture times) will permit prediction of lifetime at low loads where the 
rupture times would be very much longer. 
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(Rupture Time) 




Cyclic Loading 


Figure J1 shows some results for glass/polymer composites under cyclic 
loading. The cycles to failure have been labeled as t to simplify com- 
parison with Equations [1] and [2]. Because temperature data were not 
available, only Zhurkov's Equation [1] can be considered. The appli- 
cability of the equation is clearly demonstrated. 


Figure 12 shows cyclic loading data for metal matrix composites. Again 
the applicability of Zhurkov's Equation [1] is apparent. 


In Reference 19 considerable data were taken for unidirectional and 
angle ply graphite/epoxy composites under cyclic loading. Figure 13 
compares the results for 30-degree unidirectonal and ±30-degree angle 
plies. Similar results were obtained for the other plies of Reference 
19. The existence of parallel slopes for log r vs a has been discussed 
in the detailed development of the Zhurkov method and its modifications 
given in the next section. 


Figures 14 and .15 show the results for a and a as a function of ply 
angle for unidirectional and angle plies, respectively. Simple alge- 
braic manipulation shows the Zhurkov's y = f(a), whereas for unfilled 
polymers with parallel slopes this coefficient is zero. Here then, y 
must be associated with the fiber and/or fiber matrix interface. 


Again the utility of Zhurkov's Equation [1] is apparent. The use of 
Equation [2] is enhanced because data at only one other temperature are 
required for predictions at any other desired temperature. 


These results show that Zhurkov's approach can be used to devise ac- 
celerated tests of composites with a minimum of test data for such para- 
meters as static and dynamic tension and cyclic loading. Zhurkov's 
Equation [1] can be used at relatively high loads (short failure times) 
to obtain failure times at lads at any temperature of interest. By 
taking a few data points at one or two other temperatures, the spectrum 
of failure times can be expanded to temperatures not easily accessible. 
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D. 


THERMODYNAMICS OF CREEP RUPTURE 


The discussion in Section V.B presented the well-known Zhurkov relation- 
ship for creep rupture. When experimental results of log r(time to 
rupture) vs stress a at several temperatures are extrapolated to large 
stress, they intersect at a common pole r 0 at a critical stress o^. 
Well-behaved Zhurkov materials have poles near 10" 13 seconds. Other 
materials show poles with rg» 1Q~ 13 seconds. When these materials 
are studied in vacuum, a different behavior is observed. The earlier 
discussion showed that an entropy term, not present in well-behaved 
materials, could be justified by manipulation of the results obtained in 
air. The present discussion is a more rigorous development of these re- 
sults and shows that a premature (and probably erroneous) conclusion was 
drawn in the earlier discussion for results in vacuum. 


The following abbreviated development is after Reif.i 20 ) For a macro- 
scopic system characterized by one or more external parameters, the num- 
ber of quantum states in the energy interval between E and E + $E de- 
pends not only on the energy but also on the particular values of the 
external parameters. For simplicity, for one external parameter 
n a a (E,x). The energy of each state E n depends on the value 
assumed by the external parameter; therefore 


C4] 


For the mean value of X 


Xdx 


_3E n 

" Jr dx « dW , 


or the macroscopic work done on the system when the external parameter 
is changed quasi-statically (e.g., thermodynamical ?y reversible), is 
equal to the mean increase in the energy of the syv-«iu It can be seen 
that X is a generalized force, conjugate to the external parameter and, 
in general, may have any dimension (i.e., units). 


The change in the number of states n as a function of the work done on 
the system is given by 


( dlnfl ) _ 
9 x e 
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[53 


With the Boltzmann-defined entropy, S - k lnfy this becomes 

(iS.) *. 

ax T 

which is the desired development for the present discussion. 


Zhurkov developed his well-known creep rupture relationship from the 
time dependence of strength under unidirectional tension. He success- 
fully applied It to a variety of metals and nonmetals. The relation- 
ship is defined by 

[6] t « r 0 e(Uo - vcr)/kT 

where the time to rupture t is a function of the activation energy for 
the process U 0 less the work done ya on the system. The constant y 
depends on the nature of the structure of the material. When experi- 
mental results of log t vs stress a at several temperatures are extra- 
polated to large a, they intersect at a common point, the pole 
and the critical stress o^. A simple manipulation of Equation [6j 
shows U 0 » yct|<« The is considered to be 1/v where v is a 
vibrational frequency of bonded atoms in the solid. 

Work done at constant volume becomes a definite thermodynamic property 
of the system (i.e., a definite integral) so the work dW (Equation [4]) 

U aa #>*v»v*** 

dW a Ydcr + odY » 3(dx. 

With v being constant, Equation [5] can be written as 

[7] TaS a y(a k - a). 

From energy conservation, the heat absorbed <fQ by a system can be used 
to increase energy a E and do work. Therefore 

TaS = dQ « aE + work. 

In this quasi-static development, the work may be taken as the negative 
of the work being done on the system and, from Equation [7] 

[8] TaS * y(aj< - o) = dQ = aE - ya. 


This is just the Helmholtz "work function" relationship aA & aE - TaS 
when aA is taken as ya. With the Zhurkov-defined U Q = Y^k* su * 3 " 
stitution of o e o or a = a|/ shows U« = aE. In a given isothermal 
process, aA is the maximum {reversible) work available. From Equation 
[7], when a ~ a^ the process is adiabatic or the rupture of the speci- 
men takes place before any absorption of heat from the reservoir. When 
there is no significant PaV work, such as with solids, the Gibbs "free 
energy 1 ' aF » aH - TaS may be used for the Helmholtz relationship with- 
out loss of rigor. The association of U 0 with aH was shown earlier by 
a less rigorous development. 
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Thermodynamically Zhurkov's relationship can be written 
r a t qG^O-Y 0 ) /^T K T - ygY(a^“Cf) /kT ^^(AH-AF/kT a 

Thus at any stress level less than a|<, the time to rupture is deter 
mined by the entropy change in the solid. 


For some materials the pole T 0 is many orders of magnitude removed 
from T 0 ~ lCT 13 seconds as found for well-behaved Zhurkov materials. 
Celluloid has a T 0 ~ 10“ 4 and Lucite and Delrin have t 0 of 5 x 
10" 3 and 2.5 x 10~ 3 seconds, respectively.^ 11 * 12 ) 


When the usual Zhurkov determination for y is carried out, y is not in- 
dependent of temperature but rather is found to be y = -aj + _b where _a 
and j) are constants. For Lucite and Delrin, y becomes a zero“ at 350 and 
380K, not unreasonable values for their glass temperatures (no infor- 
mation was found for T g of Celluloid). 


For these kinds of materials, the entropy becomes 
W AS - -a(o k - 0 ) + 0) . 


Because the second term on the RHS has already been defined 
quasi-static process, the as of Equation [9] is the entropy 
a nonquasi-static process. For a quasi-static process, the 
given by 


AS = S w 


-v/ 

qs o 


dg. 

T 


as for the 
change for 
entropy is 


and for a nonquasi-static (thermodynamically irreversible) process 

as = s f - s o > f -2- 

nqs’e) 

where T now refers to the temperature of the reservoir. If a quasi- 
static path is constructed for the nonquasi-static process, the entropy 
change can be determined as 



Then Equation [9] is 
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and the extra term in Equation [9] is the discrepancy between the quasi- 
static and nonquasi-static processes. Clearly for these kinds of ma- 
terials, the rupture process is fundamentally different from well- 
behaved Zhurkov materials, with the rupture process being inherently 
nonquasi-static. 


For these materials the Zhurkov relationship can be written as 
t » t Qe“ASj/kgAS/k - i' 0 e A S/k s= TQe(^-AF)/kT _ -Y a )/kT. 


A similar result was found in a previous section. There it was stated 
that the experimentally determined t 0 » i0“ i3 seconds could be cor- 
rected by a fixed entropy term. That this cannot be strictly true is 
shown in the present development where it is found that the correction 
term ASj is a function of a. However the a of the correction term is 
quite small (0.0134 for Celluloid, 0.074 for Lucite and 0.037 for 
Delrin) so the earlier contention of a fixed term is not entirely in- 
correct. 


When Lucite and Delrin are studied under vacuum UM 2 ) the slopes of 
log t vs a differ considerably from those obtained in air. It was ar- 
gued that in vacuum the slopes at various temperature became parallel, 
i.e., no common 0 and a^. The development here shows that this is 
probably not true, and that the premature conclusion was based on a too 
narrow range of temperature used in the vacuum studies. From the de- 
velopment here it must be concluded that these materials revert to well- 
behaved Zhurkov-type materials, i.e., the finite aS^ found in air be- 
comes zero in vacuum. An obvious consideration is that this is a sur- 
face phenomenon associated with a difference in crack behavior in the 
presence of air. 


VI. CONCLUSIONS AND RECOMMENDATIONS 


The Hercules graphite/epoxy specimens under in-plane shear showed an in- 
crease in modulus after a year under thermal-vacuum exposure. Although 
the changes in ultimate strength for the AS4/3501-6 specimens are not 
considered as significant, the changes for HMS/3501-6 were large enough 
to be considered significant. When the resin content of both types of 
material is considered, the decrease in ultimate strength may be as- 
sociated with the smaller amount of resin. We recommend that this sug- 
gestion be validated by more testing with in-plane shear studies. 


In all testing of composites, a much larger number of specimens must be 
tested than has been possible to date. 


The discussion on the effects of irradiation in Section V.A shows the 
importance of the MSFC/Martin Marietta approach to such studies; namely, 
specimens must not be exposed to air after irradiation until testing is 
complete. With the paucity of literature data,U/ it is of prime im- 
portance that charged particle irradiation be subjected to all candidate 
materials for space applications. 


Zhurkov's approach to the strength of solids can be used for predictive 
modeling or accelerated testing. (8,11,12) The fundamentals of the 
method (not previously discussed in the literature) are discussed in 
Sections V.B and V.D. Section V.C demonstrate the application to com- 
posite materials. It would be of considerable value to study the sur- 
face effects noted between ambient and vacuum tests as discussed in V.B 
and V.D. Such testing would clarify the apparent discrepancy of results 
obtained under these two conditions. This would also contribute to the 
understanding of possible surface or interface effects in composites as 
discussed in Section V.C. Because the Zhurkov method is the simplest 
approach to prediction, the information developed would lead to a very 
high degree of confidence in the method. We recommend that such studies 
be considered. 


The results in Section V.C are primarily of cyclic testing analyzed by 
the Zhurkov approach. This analysis has not been previously shown. 
Cyclic tests at several temperatures would be of great importance to es- 
tablish the validity of the Zhurkov method for cyclic testing. 


For consideration at a later date it would be of value to study the ef- 
fects of irradiation on specimens while under load as discussed in Sec- 
tion V.A. 
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